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KA.TIOmL ADVISORY COMMIO?TEE FOR AERONAUTICS 

TECHNICAL NOTE 2562 

NIMERICAL DETERMINATION OF INDICIAL 
LIFT OF A TfrfO-DIMENSIONAL SINKING AIRFOIL AT SUBSONIC 
MACH NUMBERS FROM OSCLLLATORY LIFT COEFFICIENTS WITH 
CALCULATIONS FOR MACH NUMBER 0,? 

By Bernard Mazelsky 

SUMMARY 


The reciprocal equations for relating the incon 5 >re 8 si'ble circula- 
tory indlcial lift to the lift due to harmonic oscillations have been 
modified to include the noncirculatory lift associated with apparent- 
mass effects. Althou^ the apparent-mass effects are ingjulsive in 
nature in incon^jressihle flow, the lift due to apparent-mass effects in 
con^jressible flow is a time-dependent function. The corresponding 
reciprocal equatiohs for the total ccm^jressible lift are given. By use 
of the reciprocal equations for con5>ressible flow, the indlcial lift and 
moment functions due to an airfoil* s experiencing a sudden acquisition 
of vertical velocity are determined numerically for Mach number 0. 7. 

Lack of sufficient flutter coefficients prevents the calculation of 
these functions at other Mach numbers. 

Although the indlcial lift emd moment functions due to penetration 
of a shaip-edge gust may be obtained from the oscillatory tab or aileron 
coefficients by a similar analysis, sufficient coefficients are not 
available at the present. However, an approximate method is shown for 
determining a portion of this unsteady- lift function. 

When a conparison is made of the indlcial lift functions at Mach 
numbers 0.0 and 0.7^ it is noted that the growth of lift to the steady 
state appears to be less rapid for the compressible case than for the 
inconpressible case. Consequently, the calculation of the gust load 
factor at high subsonic Mach numbers utilizing the two-dimensional 
inconpressible indlcial lift functions and an over-all correction for 
conpressiblllty such as the Prandtl-Glauert factor mi^t be conservative.- 


INTRODUCTION 


A knowledge of transient flows is lnportant in many aeronautical 
problems. In the study of transient flows, two types of airfoil motions 
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have had. special sigalficance - a harmonically oscillating airfoil and 
an airfoil e2q)erienclng a sudden change in angle of attack. The lift 
function for an airfoil ejqperlenclng a sudden change in angle of attack 
and the lift function associated vith the growth of lift on an airfoil 
due to penetration of a shaip-edge gust are conmonly referred to as 
Indicial lift functions. The present paper is concerned with the use 
of reciprocal relations for detennlning for conqjressihle flow the 
indicial lift functions directly from the lift data that are available 
for the airfoil oscillating harmonically. 

The indicial lift functions have been derived for two-dimensional 
incon5>ressible flow; the function for a sudden change in angle of 
attack was derived by Wagner, see reference 1 , while the penetration 
function was derived by Kussner, reference 2 . An account of the rela- 
tions that exist between these Indicial lift functions and the lift 
coefficients for a two-dimensional oscillating airfoil is given by 
Garrick in reference 3 - 

In a recent paper by Lomax, Heaslet, and Sluder (reference 4 ), a 
method for determ in ing the indicial lift and moment function is given. 
IThile the beginning portion of the indicial functions can be calculated 
readily and the final value is considered to be the steady-state lift 
value given by the Prandtl-Glaueii; factor, the intermediate or transi- 
tion values of the indicial lift functions are difficult to obtain by 
this method and consequently numerical results are given only for a Mach 
nvcmber of 0.8 in reference 4 . 

The situation in regard to the subsonic con^jressible-flow coef- 
ficients for an oscillating airfoil is much better. Posslo (refer- 
ence 5) has fonnulated the problem in a linearized form for determining 
the lift and moment for the oscillatory case. Dietze, Schade, and 
Frazer and Skan (references 6, "J, and 8) present the lift coefficients 
of a harmonically oscillating airfoil at various Mach numbers up to 0.8 
and for various values of reduced frequency up to about 2.5. It ‘was 
felt desirable to see whether these available flutter data could be used 
in conjunction with the reciprocal relations to obtain the con5>lete 
indicial lift functions for con^jressible flow since, for dynamic- loeid 
studies, knowledge of the indicial lift functions is needed over a 
larger range of chord lengths than that given in reference 4 for Mach 
nimabers other than 0 . 8 . This paper discusses the use of the reciprocal 
relations for the case of compressible flow and presents an evaluation of 
the indicial functions for M = 0.7 for an airfoil suddenly acquiring a 
vertical velocity, for which case it was found that siifficient oscillatory 
lift data were available. 

The Indicial lift function due to penetration of a shaiq>-edge gust 
couid also be detennined by superposition provided the oscillatory lift 
coefficients were known for a -vrlde range of flap to chord ratios as well 
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as for a range of reduced freguencies. Unfortunately^ sufficient data 
are not available at the present time; however, in order to provide 
some insist as to the growth of lift for the gust case, an alternate 
method is shown for approximating most of this function. 


SmBOLS 


s, s^, 
cu 


a 


V 


c 


distance traveled, half-chords 
angular frequency 
forward velocity of airfoil 
chord 


k 

L(s) 

M(s) 

P 

h 

0 •• 

h, h 
M 

6(s) 


reduced- frequency parameter 



lift per unit length of span 

moment per unit length of span about quarter- chord point 
density 

anqjlitude of vertical displacement, half-chords 
first and second derivatives of h with respect to s 
Mach number 

inq)ulse function, B(s) = “> for s = 0, B(s) = 0 for s ^ 0, 

poo 

and / 6(s)ds=l 


k, (s) indicial lift function for an airfoil experiencing a sudden 

acquisition in vertical velocity as used in equation 

L = -JtpcV^hk]_(s) 

U velocity of sharp-edge gust 

k^Cs) indicial lift function due to penetration of a sharp-edge 

gust as used in equation L = -itpcVUk 2 (s) 
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k 

C(k) 

Cc(k) 

M(k) 

f(k) = 
m(k) = 
Z2 

M2 

A, B, C 
Matrix 

L ] 

{ } 


indlcial moment function for an airfoil experiencing a 
sudden acquisition of vertical velocity -vdiere moment is 
taken atout quarter-chord position 

center-of-pressure location in percent chord from leading 
edge for an airfoil experiencing a sudden acquisition of 
vertical velocity 

Theodorsen’s circulatory lift function (reference 9), 

F(k) + iG(k), as used in equation L = -jtpcV^he^® jj-ikC(k^ 

conplex conpressihle-flow oscillatory lift coefficient -which 
includes hoth circulatory- and noncirculatory-lift com- 
ponents, Cc(k) = lised in equation 

L = -itpcV^e^hj^ikCc(k^ 

in-phase conqwnent of the conplex moment corresponding to 
■t±ie F^(k) lift coefficient 

F^(k) - F^(~) 

M(k) - M(«) 

in-phase and out-of-phase lift coefficients associated -with 
translation of the airfoil as used in equation 

L = itpcV^e^® + iZ2) 

in-phase and out-of-phase moment coefficient about quarter- 
chord position associated -with translation of the airfoil 

as used in equation M = jtpc%^^® ^(^1 


! cons-tanta 
no-tatlon: 


rectangular matrix 


column matrix 
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ANALYSIS AND EESULTS 

Superposition Integrals for. the In;coii5)ressible 
and Con5)ressible Cases 


The reciprocal eq.uations idiich are essentially a special form of 
the supeip)osition integral have been derived for the incompressible 
circulatory lift by Garrick (reference 3)* This section first discusses 
an extension of Gairlck’s work by incoip>orating the potential lift due 
to apparent-mass effects into the reciprocal equations for the incom- 
pressible case. ' The reciprocal eq,uations for the con5)ressible case are 
then indicated and the evaluation of the Indlcial lift and moment for 
an airfoil experiencing a sudden acquisition of vertical velocity is 
determined for the case of M. = 0.7. 

Incompressible case .- In. the incon5)resslble case, the lift in 
unsteady motion has been separated into two parts - a noncirculatory 
lift Tdiich is commonly referred to as an apparent-mass effect and a 
lift due to the circulation about the airfoil. The reciprocal eq.ua- 
tions for incon5)ressible flow are expressed only in terms of Theodorsen's 
circulatory- lift function C(k) and Wagner’s indlcial- lift function k]^(s) 
and can be expressed in two separate forms as follows: 


(s'>0) (1») 

and 

k^(s) = 1 * f r di (= > 0) (lb) 


where F(k) and G(k) are, respectively, the in-phase and out-of-phase 
con^jonents of circulatory lift .on a harmonically oscillating airfoil defined 
by the equation 


L = -«pcv\e^[llkC(k2 


where 


C(k) = F(k) + iG(k) 
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Tte kQ^(s) function Is defined such, that the lift on an airfoil experi- 
encing a sudden acquisition of "vertical velocity may he expressed hy "the 
following equation: 


L = -itpcV^fik2^(s) (2) 


Two relations that can he oh"tained from equation (la) are (l) that 
the "value of F(k) at k = 0 is equal to "the asynptotic "value 
of k]^(s) as s tends to infini"ty, -these "values corresponding to "the 
steady-state lift, and ( 2 ) "that the initial "value of k2(s) is equal 

in magnitude to the asynptotic "value of F(k) as k tends to infinity. 
This second relation is not readily e"vident from equation (la); con- 
sequently, a ma"thematical proof of this second relation is shown in 
appendix A. As a result of these relationships, the initial "value and 
the asynpto"te of "the k]_(s) function can he determined from the given 

F(k) function. A knowledge of "these end points is valuable in a 
numerical solution of the reciprocal equations. 

As noted pre"viously, equations (la) and (ih) express only the 
circulatory part of the lift and, therefore, do not account for the 
apparent-mass effects. An expression for the apparent-mass effects for 
the indlcial case, however, may he obtained as follows. 

The lift due to apparent-mass effects alone for any arbitrary 
vertical motion of "the airfoil when written in terms of the nondi- 
mensional displacement h is given in reference 9 s-s 




( 3 ) 


When this equation is applied to "the case of an airfoil experiencing a 

sudden vertical velocity, h(s) becomes inpulsive in character. The 
magnitude of this inpulse may he defined hy 


h(s) = h 6 (s) 


where h is "the lns"tan"taneous vertical velocity acquired hy the airfoil 
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and the function 6(s) Is .defined such that 6(s) = 00 ' for s = 0, 
5(s) = 0 for B f 0, €ind 


/ CO 

5(s)ds = 1 

■ 00 


• • • 

When h8(s) is substituted for the iB 5 >ulsive sinking acceleration h(s) 
in eguatlon (3), the resulting expression for the lift due to apparent 
mass alone becomes 


-■iiipulse 


= -jtpcV^h 


S(s) 

2 




If the inpulsive part of the unsteady-lift functions is designated by 
ki (s)^TTipni pp aiid if the form of equation (2) is retained, then the 

expression for the inpulsive lift may be -written 


^inpulse ® ••^PcV^hkp(s) (5) 

The value of kjL(s) due to apparent-mass effects in the incanpresslble 
case may then be determined by a conparlson of equations (4) and (5) to 
give 


' 5 ( 6 ) 

Addition of -the k]^(s) function due to circulation, equations (l)., and 
the l^i(s)fji^u 2 gQ function, equation (6), then gives the expressions 
for the lift for -the Indlcial case as follows: 

l^l( = )tctal - f/J d:. . 1 S(b) 


(s ^ 0) (7a) 
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and 


l^l(«)total = 1 G(k) cos _ ks , at + 1 5(g) (g > q) (7b) 


CQfflpresslTjle case .- While In the incongpressihle case it has heen 
found convenient and natural to consider the circulatory and non- 
circulatory lifts separately (this separation is especially desirable 
because of the in^julsive nature of the noncirculatory part), these 
circumstances do not exist in the con^iresslble case. The perturbation 
velocities are finite in the con^resslble’ case (since speed of sound is 
considered finite) and, consequently, the effects of the pressure dis- 
turbances in the "whole flow field are not felt instantaneously as for 
incongoresslble flow. Consequently, the lift due to apparent-mass 
effects would not be In^julsive but would remain finite and time dependent. 
Separation of the lifts is not necessary in the compressible case and, 
moreover, is not convenient since the available results for the com- 
pressible case are presented numerically only for the total lift. The 
reciprocal equations for the total conqpressible lift are, therefore, 
much simpler than the equations for the incon5)ressible case since the 
inpulse function may be disregarded. 

For con^jresslble flow, the case is exactly analogous to the treat- 
ment of indicial admittance given by Bush in reference 10. In this 
reference, equations are given which relate the indicial admittance to 
the in-phase and out-of-phase responses of a system subjected to an 
oscillatory forcing function. In notation similar to that- used for 
incoii5)3ressible flow, iiese ejqpressions may be written as follows 

^l(s) total = f k ^ (b > 0) (8a) 

and 

° m (b> 0) (8b) 

where F^,(k) and G(.(k) are, respectively, the in-phase eind out-of- 
phase lift coniponents on an oscillating airfoil end are defined by the 
equation 
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L = -itpcV%ie^® -il£Cc(k^ 


(9) 


■where 

C^(k) = F^(k) + iGc(k) 


Equations (8a) and (8b) are appUcal^le for the con5>ressible case pro- 
viding -the function Cj.(k) is continuous and finite in -the interval 0 

to 00, From the physical conditions for the conpresslble case, it can 
reasonably be assumed that these functions adhere to those conditions. 

It is of interest to note that, -vdien these equations are applied to 
lncon5)ressible flow, they apply only to -the circulatory part of the lift 
as given by equations (la) and (lb). In this case, "the factor Fj.(0) 

becomes, by definition, sluqjly unity. For conpressible flow, the fac- 
tor F^(0) acquires a value given by the Prandtl-Glauert factor 5= . 

l/TT^ 

The existing coefficients for the translatory case -vdiich are to be 
used in the reciprocal relations have been given in -various forms begin- 
ning with Possio’s -work. Ihe form in -which all the results are commonly ■ 
put is 


L = jtpcV^e^® + IZg) (10a) 

and 

M = itpc^^e^ |(M3^ + IM2) (lOb) 


where M is the moment acting about the quarter-chord point. In order 
to make convenient use of the reciprocal relations given by equations (8a) 
and (8b), it is first necessary to convert the expression for lift given 
by equation (lOa) to a form similar to the one used in equation (9). When 
the expressions (9) and (lOa) are con5>ared, the following correspondence 
is seen to exist for conpresslble flow; 


Fc(k) 


Z2(^) 

2k 




(11a) 
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and 


Gc(k) 


2k 


(lit) 


Summary of Available Flutter Coefficients and Related Data 
for Subsonic Compressible Flow 

The two-dimensional conjiresslble flutter coefficients for the real 
and imaginary parts of the lift and moment oscillatory coefficients for 
sinking motion are given in table I for three Mach numbers, M = 0.5j 
M = 0.6, and M = 0.7. This table represents a sumnary of the results 
of three authors where the flutter coefficients by each have been con- 
verted to the Z]_, Z 2 , M]^, and M 2 form for use in equations (lOa) 

and (lOb). The values of and Z 2 correspond to the lift coef- 

ficients indicated in equation (lOa), and the values of M-j^ and M 2 

indicated in equation (lOb) correspond to the moment taken about the 
quarter-chord position. The values of these coefficients were obtained 
from three sources, eind the range of reduced frequencies taken from each 
source is indicated in table I. Data for large values of reduced fre- 
quency do not appear to exist. The accuracy of these data for reduced 
frequencies lower than approximately 1 is better than the accuracy of 
the data at hi^er values of reduced frequency (see reference 8), but a 
reasonably good solution for the indicial lift and moment functions is 
stlid. possible since this decrease in accuracy at the higher reduced 
frequencies will be shown to' have little effect on the determination of 
either the unsteady lift or unsteady moment functions provided the value 
of the flutter coefficients at ilie infinite reduced frequency is known. 
Fortunately these values can be determined from the initial value of the 
indicial functions (see appendix A). If the flutter coefficients are 
known up to a reduced frequency where they have practically the same 
value as the flutter coefficient at the infinite reduced frequency, then 
they are known for a sufficient range of reduced frequencies for deter- 
mining a reliable solution of the indicial functions. This condition 
exists for M = 0.7 since values of the flutter coefficients at a 
reduced frequency of 2.5 have been found to agree quite well with the 
values at the infinite reduced frequency. However, values of the flutter 
coefficients at M = 0.5 and M = 0.6 are known only for reduced fre- 
quencies up to 1. Also, the flutter coefficients at these Mach numbers 
must be known for hi^er values of reduced frequency than for M = 0. 7 
since the flutter coefficients approach the values at the Infinite fre- 
quency at a much slower rate. Consequently, the only Mach number for 
which sufficient data were found available for obtaining the indicial 
functions is M = 0.7. 
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In order to aid in the numerical solution of the reciprocal equa- 
tions^ the end points of the flmction should he determined. 

These end points may he detennined independently of the flutter coef- 
ficients. The value of 'vdiich represents the steady-state lift 

can he given hy the Piandtl-Glauert factor When an atten^jt is 

\/l - 

made to calculate this end point from the and Z2 coefficients 

given in table I, an indeterminate form is found; however, the magnitude 
of this indeterminate quantity 'vdien properly evaluated agrees very well 

with the factor — . The value of F(,(n>) -which corresponds to the 

]/l - 

value of k-^Cs) a‘t s = 0 (see appendix A) can he determined hy the 
following equation 


Numerical Solution of the Reciprocal Equation 

The indlcial lift function k3^(s) can he determined from either of 
equations ( 8 a)‘ and ( 8 h) -where the coefficients Fg(k) and Gg(k) are 

determined from the oscillatory coefficients hy equations (lla) and 
(Uh). It has heen noted thatj of the data presented in table I, the 
only Mach number for -which sufficient data are available for providing 
a reasonable estimate of the indlcial lift and moment on the airfoil ' 
is M = 0.7 since values of the flutter coefficients are not given for 
values of k > 1 at the other Mach numbers. A plot of the Fp(k) and 

Gj,(k) functions for this Mach number is shown in figure 1 . Althou^ a 
smooth curve can he dra-wn for the F,,(k) function at the hl^er reduced 
frequencies, -the data for the G^,(k) function appear to he erratic in 

this region. Consequently, the latter function is represented hy a 
broken line a-t the hi gher reduced frequencies. Also, the numerical 
e-valuation of the k3_(s) function hy equation ( 8 h) is more difficult 

than hy equation ( 8 a) because -the integrand in equation ( 8 h) is inde- ' 
terminate at k = 0 . In order to evaluate -this indeterminate form, the 
value of the derl-vative of G^,(k) at k = 0 must he known. Since the 

function G(,(k) is not given in a closed form, this derl-vative would 

he difficult to determine accurately in the conq)resslhle case. Because 
of the erratic nature of the G^.(k) function at the hi^er -values of 

reduced frequency and because of the presence of this indeterminate form 
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of k = 0, equation (8a) is used tlirou^out the ensuing 

analysis for determining the indicial functions. 

The graphical solution of equation (8a) may he sin^jlified consider- 
ably hy incorporating in the solution of this equation the known value 
of the flutter coefficient at the infinite reduced frequency hy making 
the following suhstitution for the flutter coefficient Fc(k): 


or 


f(k) = P^.(k) - F^(«.) 


(13a) 


F^(k) = f(k) + F^(c») 


( 1313 ) 


The suhstitution of equation (131>) into equation (lla) leads to an 
alternate equation for kQ^(s): 


ki(s) = F^(oo) + SM. 


sin ks 


dk 


(11^) 


The graphical solution of k]_(s) is thus simplified since the integral 

in this expression can he evaluated more readily than the integral in 
equation (8a) because the Integrand approaches zero much more rapidly. 

The F^(k) function for M = 0.7 is calculated from the flutter 
coefficient Z 2 by equation (lla) and is shown in figure 1 together 

with the value of the asymptote calculated hy equation (12). The trans- 
formed function f(k) as obtained hy equation ( 13 &) is shown in fig- 
ure 2. A plot of the Integrand in equation (l4) for several values of s 
is shown in figure 3- These Integrands were Integrated hy means of a 
planlmeter. The k^(s) function thus found hy means of equation (l4) is 

shown in figure 4 along with the k]^(s)- function for M = 0.0. In this 

same figure a part of an Independent solution for M = 0. 7 given in 
reference 4 is also shown. 

It will he useful to fit the k 2 (s) function for M = 0.7 given 

in figure 4 hy some analytical function. Since the exponential function 
has a simple operational equivalent and has been found convenient in 
approximating the k-]_(8) function at M = 0.0 (see reference U), a 
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limited series of such functions was chosen to approximate the 

k]_(s) function at M = 0.7. The function found to fit this curve quite 

■well is 


kp(s) = 1.4(1 - + 0.4l9e"°*^^®) (15) 


and is also shown in figure 4. The conre spending approximate expressions 
for the harmonically oscillating airfoil can' he found from eq.uatian (15) 
in a similar manner shown in reference 11. The expressions for 
and G(j(k) are 


F.(k) = 1.4 


1 - 


0. ^64k'' 


0 . 40 ^^ 


0.4l9k^ 


2 2 2 2 2 2 
(0.0536) + k (0.357) + k (0.902) + k 


(16a) 


and 


G_(k) =1.4 


0.01951k 


0.l446k 


Siz.m9X 


2 2 2 2 ' 22 
(0.0536) +k (0.357) +F (0.902) +k 


(l6h) 


Both of these expressions are plotted in figure 1 for conparison with 
the data obtained from the flutter coefficients. Note that the approxi- 
mate expression for the F^(k) function is in fairly good agreement for 
all of the k values considered. However, relatively poor agreement 
exists for the G^(k) function at the hi^er values of k. This com- 
parison is to he expected since, as was previously cited, the Gj,(k) func- 
tion obtained from the flutter coefficients appears to "be q,uite erratic 
at the hi^er values of k. 

The solution of the reciprocal eq.uation has thus far been shown for 
the calculation of the k-j^(s) function. Since the reciprocal equations 

are also applicable to the determination of the imsteady moment n^(s) 
due to a sudden change in airfoil vertical velocity, the m^(s) function 
taken about the quarter- chord point can be shown to be given by 
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f 2 r” M(k) sin ks ,, 

■ Ijo -^-"4 — 


( 17 ) 


■where M(k) is considered, herein to he determined from the flutter coef- 
ficient M 2 hy the following equation: 


M(k) 


Mg 

'K 


( 18 ) 


A plot of the function M(k) for M = 0.7 is shown in figure 5 together 
•with the -value of M(<») determined hy 


M(“) 


1 

2idl 


( 19 ) 


In a manner similar to that shown for the unsteady- lift case, the func- 
tion M(k) is then transformed in order that equation (IT) can he 
e-valuated graphically hy the following substitution: 


m(k) = M(k) - M(“) (20a) 

or 

M(k) = m(k) + M(oo) (20h) 


4 

A plot of the function m(k) for M = 0. T .is shown in figure 6. The 
substitution of equation (20b) into equation ( 17 ) leads to the expression 

m^(s) = M(») + f J di (21) 

Plots of the integrand in equation (2l) are shown in figure 7 for several 
values of s, and -the n^(s) function e-valuated hy this equation is shown 

in figure 8. 

If the indicial lift function k^(s) ' and the indicial moment func- 
tion n^(s) are known, the variation of indicial center-of -pressure 
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location In percent chord can be calculated by means of -the following 
equation: 



m-i (s)lOO 

= 25 - 

kj^(s) 


( 22 ) 


A plot of the center- of -pressure location is given in figure 9 and, for 
con5>arison, the case for M = 0.0- is also shewn. 


Approximate Considerations for the Indiclal Lift hue to 
Penetration of a Shaip-Edge Gust 

Because of the limited number of flutter coefficients available for 
flap motions, as noted in the introduction, the determination of the 
indicial lift function for penetration k2(s) by use of the reciprocal 

equations is not possible at the present time. An alternate and perhaps 
semlratlonal method for determining a part of the penetration func- 
tion k2(s) from the kp(s) function on the basis of the relatians 
between these two functions for the inconpresslble case is suggested as 
follows. The relations between the k^^Cs) and k2(s) functions for the 

lncoD5)ressible case are given in reference 3 and are as follows: 



ds]^ + i^s(2 - s) (0 < s < 2) (23a) 



(s >2) (23b) 



The first term in equation (23a) is associated with only the circulatory 
part of the k]^(s). function -vdiile the second term of the equation arises 

from the inpulslve term in the k]_(s) function, and hence is associated 

with the apparent-mass effects. This second term is known not to apply 
for conpressible flow, since for the k2^(s) function the apparent-mass 

effect is no longer lu^ulsive in character but rather is time dependent. 
Therefore this equation cannot be applied to the compressible case. With 
regard to the applicability of equation (23b), it is assumed that this 
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equation can evaluate tlie circulatory part of the lift for the k2(s) func- 
tion provided that only the circulatory part of the k]_(s) function is 
used. Althougjh the kj_(s) function calculated herein for compressihle 

flov contains hoth the circulatory and noncirculatory lifts, the following 
argument may he used as a hasls for the assumption that for values of s 
greater than approximately h, the function can he considered, for 

practical purposes, to he due to circulation only. If the assunption is 
made that the lift due to apparent-mass effects and the lift due to circu- 


lation act always at chord position and chord position, respectively, 
the moment function n^(s) shown in figure 8 must arise from the apparent- 


mass effects only. This figure shows that for practical purposes the 
moment function, referred to the quarter-chord position, has decayed to a 
negligible value for s greater than approximtely 4.' Thus the portion 
of the k-j_(s) function in figure 4 for s greater than approximately 4 


may he associated irf-th the lift due to circulation only. Inspection of 
equation (23h) shows that, if the evaluation of the k2(s) function is 
confined to values of s greater than 6, only the portion of the 
k^(s) funqtion for s greater than approximately 4 (the portion attributed 

to circulation) will he used. A numerical solution of equations (23a) 
and (23b) is given in appendix B. The k2(s) function for M = 0.0 was 

calculated by using the k2(s) function for M = 0,0 shown in figure 4 


as a means for estimating the accuracy of a numerical solution given by 
equation (B9) in appendix B. The values obtained were within 1.5 percent 
of the known Kussner function. For the conpressible case the k2(s) func- 
tion for s greater than approximately 6 was evaluated by utilizing the 
kj^(s) function for M = 0.7- The results of this calculation are shown in 

figure 10. Althou^ the solution for the conpressible case is approximate, 
it is Interesting to note that for s > 6, "vdiere the solution was reasoned 
to be applicable, the grovrth of lift is less rapid for the conpressible case 
than for the inconpressible case. This phenomenon would be expected since 
the growth of lift due to the kQ^(s) function, which was used to determine 

the k2(s) function, shows this same effect as evidenced by figure 4. 


POSSIBLE EFFECTS ON GUST LOAD FACTOR 


The fact that the growth of lift for the k]_(s) and k2(s) functions 

is less for the M = 0.7 case than for the M = 0.0 case may have some 
significance in the determination of the response of an airplane to gusts. 

The current design procedure for the loads due to gusts on an alzplane 
utilizes the two-dimensional inconpressible- flow k]^(s) and k2(s) functions 
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and an over-all approximate correction for con 5 )ressi‘billty. If in ttie 
calculation of the gust load factor at about M = 0.7 the corresponding 
k]_(s) and kgCs) functions for M = O.7 vere used Instead of the 

Incon^jresslble k]^(s) and k 2 (s) functions with the over-all approxi- 
mate correction for coirpressibillty, then the resulting load factor 
would probably be less than that found by using the Incon^jressible 
k-j^(s) functions. This further suggests the desirability of obtaining 

the k 2 (s) function for the con 5 )lete range of s and of obtaining the 
k 2 (s) functions for other values of Mach number so that a greater 

insight may be had as to the effect of con^jressibllity on the gust load 
factor. 


CONCLUDING REMABKS 


The reciprocal equations which relate the incon^jresslble circulatory 
lift to the lift due to harmonic oscillations have been modified to 
include the noncirculatory lift due to apparent-mass effects. Nhile the 
lift due to apparent-mass effects for the indlcial case is known to be 
impulsive for the incompressible case, the lift due to apparent-mass 
effects for the corresponding compressible case was found to be finite 
and a time-dependent function. 

The reciprocal equations are evaluated numerically for determining 
the indiclal lift and moment on an airfoil experiencing a sudden acqui-^ 
sltion of vertical velocity at M = 0.7* Lack of sufficient flutter 
coefficients prevents the calculation of these functions at other Mach 
numbers . 

The growth of lift to the steeidy- state value on an airfoil experi- 
encing a sudden acquisition of vertical velocity was found to be less 
rapid for the compressible case than for the incompressible case. 

An approximate calculation is made for a portion of the indiclal- 
lift function due to penetration of a sharp-edge gust at M = 0.7. It 
is noted that the growth of lift to the steady- state value for this - 
indiclal lift function was also less rapid than for the incompressible 
case. As a consequence of these phenomena, the calculation of the gust' 
load factor at hi^-subsonic Mach numbers utilizing the two-dimensional 
incompressible indlcial lift functions and an over-all correction for 
compressibility such as the Prandtl-Glauert factor mi gbt be conservative. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., July 30 , 1951 
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APPENDIX A 


PROOF FOR lim XtCs) = F„(oo) 
s->40 

The function ^-^(8) can he expressed as a function of F^,(k) hy 
the following equation; 


, 2 r Fc(k) Bto ka 
i -hJq k 


dk 


(Al) 


provided Fj.(k) is a hounded and well-hehaved function. The value of 
k]_(s) as s tends to zero positively may by evaluated by taking the 
limit of both sides of equation (Al) as s — >+0; 


lim kp(s) 
s — 


lim 2 r°° Fc(^) sin ks 

s— >+0 ^^0 ^ 


(A2) 


Before the limits in equation (A2) are evaluated, it is convenient to 
make the following substitution for Fc(k): 


F^(k) = f(k) + F^(oo) (A3) 

Substitution of this expression for F^.(k) into equation (A2) leads to 
the following expression for k]^(s): 

Ita k,(s) = Um s.( r dk + F^(.) r dk) (Al.) 

g_^0 s— ^ ^ / 
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Althou^ the second integral in equation (Ah) can he evaluated readily 
to give 


lim / sln^ ^ ^ 1 (A5) 

s — >+ 0'-'0 ^ ^ 

the evaluation of the first integral in equation (Ak) is not readily 
evident; however it can he shown (see lemma) that 


llm r" sto ka ^ ^ Q 


since 


(A6) 


lim f(k) = 0 
k — ^ 


As a consequence of equations (A5) and (a 6) the value of the limit in 
equation (Ak) is 


lim ki(s) = Fc(«>) ’ (AT) 

s — >+0 


Lemma; 

Prove that 


lim 
s — )+0 



f(k) sin ks 
k 


dk = 0 


°(A8) 


provided f (k) is a hounded and integrahle function within the limits of 
integration and 

lim f(k) = 0 
k — ^00 
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If the folloadJig suhstitution is made for ks 


ks = X 


an alternate estpression for equation (A8) is had 


lim 
8 — )+0 



sin X 

X 


dx = 


0 


(A9) 


Since is hounded and integrahle within the limits of Integration 

and since the integral 

poo 

/ sinx ^ 

Jo X 

is absolutely convergent, then equation (A 9 ) is uniformly convergent 
for s ^ 0 within the limits of integration. Consequently, there is 
a positive number X such that 





^ dx 

Jx 

)/ X 



e 


(AlO) 


where e is an arbitrary number chosen as small as desired and the 
number X is independent of the value of s. Also a number Xq, inde- 
pendent of s, may be chosen so that 


f 


dx 


<i s 


(All) 


However, 
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■where 


Xo 1 ^ < X 


It can be shown that 


therefore 


L 


sin X 


dx 


< « (q. > P > 0 ] 


f f(f)^ a. < 


thus a -value of S cein be chosen as small as desired so that 



is%S) 


Combining the results given by equations (AlO), (All), and (A12) 
the equation 



thus, in the limit as S — >0 positively 


Urn r dx = 0 

s— >+oJo ^ 


or 


lim 

s-^+0 



f(k)sln ks dk 


= 0 


(A12) 

yields 

(A13a) 

(A13b) 


k 


■ (A13c) 
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APPENDIX B 


NUMERICAL EVALUATION OF INCCMPEESSIBLE-FLOW INDICIAL LIFT FUNCTION 
DOE TO PENETRATION OF A SHARP-ED(^ GUST 


The indicial lift function due to penetration of’a sharp-edge 
gust k2(s) expressed in tenns of the Indicial lift function k3_(s) 

is given hy equations (23a) and (231>). Althou^ a numerical evaluation 
of the first term in equation (23a) may he performed -with the aid of 
Siupson’s rule, the numerical evaluation of the integral in equation ( 23 h) 
is difficult in that the integrand is infinite at S]^ = 2 . In order to 

oveix:ome this difficulty, a paraholic-arc segment can he fitted to the 

function - ®]_)^ between the limits Sp = ^ = 2 . Because 

of the similarity between equations (23a) and (23h), the numerical solu- 
tion of both equations can be written in one eq.uation with the aid of 
matrix notation. For simplicity, the numerical solution is calculated 
at integral values of s. However, for increased accuracy the interval 
of integration is taken at quarter multiples of s. In view of these 
considerations, the numerical solution for the k2(s) function may be 

performed in three steps; 

( 1 ) Integration of integrand for the limits = 0 to s^ = 1 

With the aid of Slnpson’s integrating factors the following integral 
is evaluated, a quarter multiple of s being assumed for the interval 
of integration; 



(2) Integration of integrand for the limits s^^ = 1 to s^^ = ^ 

In a manner similar to that shown for the evaluation of the integral 
in step 1 , SLupson’s integrating factors can be utilized to calculate the 

same integral over the limits s^ = 1 to ^ as follows; 
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d.S- <« i i 

1 It 14. 


i(l)kl(s - 1) + 


k 

3 




(B 2 ) 


( 3 ) Integration of integrand for the limits ^ to s-j^ = 2 

As noted previously the integration for these limits is not possible 
hy means of the sin 5 )le Simpson’s factors since the integrand is infinite 
at = 2. With the change in variable = 2 - a, a parabolic-arc 

segment is fitted to the function - s^^ between the limits s-j^ = 

to Bj^ = 2 by use of the following equation: 


k^Cs 


2 + a) <« A + Bu + 


Ccr^ 


(B3) 


A set of s imu ltaneous equations relating the function k^Cs - 2 + cr) with 
the constants A, B, and C in the parabolic-arc segment may be obtained 
by evaluating equation (B3) at the pointa o = 0, ^ ^ ~ 2 

follows: at o = 0, 


at 



at a 


1 

2 



(Bh) 


(B5) 


(B6) 


rolOJ 
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^Jtien equations (b 4), 
tlie constants A, B, and 


(B5), and (b 6) are solved simultaneously for 
C, the resulting expressions are 


A = k3_(s - 2) 


B = -6k2^(s - 2) + Skj^l^s 
C = Skj^Cs - 2 ) - l6k;Lfs 


1 ) - - i) 

i) * - 1) 


> (BT) 


If the expression for kj_(s - 2 + o) given hy equation B(3) is substi- 
tuted into the integrand inider consideration together with the values 
A, B, and C determined hy equations (B7), the integral can he 


of 


evaluated in closed form in terms of the unsteady- lift function kj^^s - s-j_^ 


from Si 


^ 1=2 
as follows: 


to SjL = 2 . The results of this closed- form integration are 


"/s/a ^ 


3\/3 - 2k 


»(• - i) • ■>(• ■ 8 


(B8) 



D 


i 


Combining equations (Bl), (B2), and (b 8) allo-wn the solution of the two equations (23a) and (23b) 
for the ko^®) function to be written in a single matrix equation for integral multiples of s 
as follows: 


'’ke(i) 


J:- 

12 

k£(2) 

1 

1 

— K TV 


12 

ite<3) 

1 

SC 

0 

1 

k(iO 

' 

0 


Ts-WJ^ l/W 


T 

t 


'M 

13 


3 


ZL 


0 0 0 0 0 


u U V V V 




1£ 

0 


13 


£ 

21 


0 0 0 0 


0 -3^/3_ac 1 

^ 6 9 6 15 10 ZL 


kl(0) 

‘>(1) 

Mi) 

kl(l) 

M'i) 

J.\ K/ 

ki(e) 

ki(^) 

"i(^) 

kl(3) 

"l(4) 


[A 

1C 


(B9) 
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In figure 10 the solutions for the k^Cs) function are given 

for M = 0.0 n.nri M = 0.7* No discemahle difference "was present idien 
the calculated values for M = 0. 0 were con^jared with the Kilssner 
function. 
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TABLB I.- CQEFFICIEHTB FOR THE LIFT AND MCHfflO; 



Zl]_ for - 


Z2 for 


M = 0.5 H =■ 0.6 M a 0,7 M - 0.5 M - 0.6 M = 0.7 M « 0.? 


00»H7l .005231 .00711 .043851 .04699 
.02006 .08308 
.03399 .1183 
.04722 .1503 
.0589 .1796 

.0945 .3057 

.0991 .4175 

.0827 .5286 

.0538 .6447 

.0187 .7698 

-.0154 .9071 

-.0449 
-.0874 
-.39295 
-.6521 
-.6083 


1324 


1651 


1941 


3186 

4332 

-.01341 

-.02944 

5523 

6820 

-.0518 

-.0805 

8229 

-.1160 


for - 


M2 for 


Bfiference 


m 

M - 0.7 






0 

0 

0 

0 




-.00026 

.000005 

,00001 






.00004 





00158 

-.00208 

.00010 

.00020 




00274 

-.00355 

.00021 





0042 


.00034 




^ 6 

0155 

-.0186 

.00163 





0336 

-.0393 

.00389 

.00712 

.01346 



0585 

-.0672 

.0074 

.01354 

.02596 



0902 

ioi4 

.0126 

.02324 

.04482 



1260 

-.1402 

.0201 

.03752 

.07175 



1742 

-.1613 

.0308 

.0581 

.1090 



2254 

-.2178 

.0438 

.08185 

.1578 

1 

1 

3455 

-.2798 

.0881 

.1718 

.2665 


f ^ 





4896 

“\ 






7n R“i 









^ 8 


k'OL’r 

ViliWI 


1.218 
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8.- Confiariaon of the Indicial moment function due to a 

[en acquisition of vertical velocity at M = 0,0 and M ■ 0,7. 
sent taken about quarter chord, ) 
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Figure Conparison of center-of-pressiire variation from leading edge of 


X J A . 
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and M = 0.7. 
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Figure 10,- Coinparison of the indlcial 1: 


a sharp-edge gust at M = ( 







